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The  Behavior  of  Gridded  Spherical  and 
Planar  Electron  Probes  in  a 
Non-Maxwellian  Plasma 


1.  CURRENT  CALCULATIONS 


Since  the  pioneering  efforts  of  Mott-Smith  and  Langmuir,  many  investiga- 

2 

tions  have  analyzed  the  response  of  electrostatic  probes  in  stationary  plasmas. 
Because  of  their  sensitivity  to  low  energy  particles  (<1  eV),  electrostatic  probes 
of  various  geometric  configurations  have  been  used  to  measure  ionospheric  param- 
eters. In  the  coming  year,  several  experiments  aboard  the  SCATKA  satellite  will 
attempt  to  measure  the  thermal  properties  of  plasmas  at  ~6Rg  above  the  equator. 
At  such  great  altitudes,  it  is  possible  that  the  distribution  function  of  thermal 
electrons  varies  from  the  Maxwellian  distributions  commonly  found  in  the  lower 
ionosphere.  In  anticipation  of  such  a variation,  we  here  consider  the  responses 
of  gridded  spherical  and  planar  electrostatic  analyzers  to  a type  of  non-Maxwellian 
plasma  already  encountered  in  several  regions  of  space.  The  probes  that  we  con- 
sider here  are  made  up  of  either  a collecting  sphere  surrounded  by  a concentric 
wire  mesh  grid  (Figure  la),  or  a planar  collector  with  circular  gridded  aperture 
(Figure  lb).  The  collectors  are  held  at  sufficiently  high  positive  potentials  to 


(Received  for  publication  15  March  1978) 
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Figure  la.  Velocity  of  a Particle  Relative  to  a Gridded  Spherical  Electrostatic 
Analyzer  of  Radius  Rq.  (e.  and  (e,^)  are  the  colatitude  and  azimuth  angles  in 
configuration  and  velocity  spaces,  respectively 


Figure  lb.  Velocity  of  a Particle  Relative  to  the  Circular,  Gridded  Aperture  of 
a Planar  Electrostatic  Analyzer.  The  aperture  radius  is  and  (e.X)  are 

the  colatitude  and  azimuth  angles  in  configuration  and  velocity  spaces,  respectively 


screen  ambient  thermal  ions  from  detection.  The  grids  may  be  either  swept  or 
stepped  in  voltage.  Similar  instruments  have  been  flown  on  several  Air  Force 
satellites.  Equations  appropriate  for  spherical  probes  are  developed  explicitly  in 
the  text  and  those  for  planar  geometries  in  Appendix  A. 

The  current  (I)  measured  by  a gridded  spherical  electrostatic  probe  with  an 

applied  retarding  voltage  (V)  due  to  a Maxwellian  electron  gas  of  temperature  (T) 

. 3 

IS 

.aX 

I(V)  = TTR^  ^ n q c e (1) 


where:  fi  and  are  the  transmission  coefficient  and  radius  of  the  grids;  n,  q, 
and  c are  the  electron  density,  charge  and  mean  thermal  speed  (SkT/^rm)  ' , 
respectively;  and  m is  the  mass  of  an  electron  and  k is  the  Boltzmann  constant. 

T is  calculated  by  examining  the  current  characteristic  (In(-I)  versus  V curve) 
for  a linear  portion  and  using  the  relationship 

T = - ^ Id  In  (-I)/dV)"^  . (2) 


In  the  analysis  of  electrostatic  probe  measurements  taken  in  the  ionosphere, 
it  is  commonly  assumed  that  the  plasma  is  Maxwellian.  In  the  lower  reaches  of 
the  ionosphere  where  collision  times  are  short  in  comparison  to  the  measurement 
time,  the  Maxwellian  assumption  is  justified  by  the  H-theorem.  In  the  topside 
ionosphere,  the  assumption  is  justified  by  experience.  Frequently,  however, 
plasmas  found  in  deep  space  have  long  relaxation  times.  In  such  cases,  a Max- 
wellian interpretation  may  be  invalid.  For  example,  the  electrons  measured  in 
the  plasma  sheet^  and  magnetosheath ^ have  "kappa-type"  distribution  functions  of 


^K<v)  = ^ 

w 


r(K  + 1) 


1 


(jrK)^/^  r{K-l/2)  |l  + 


2 TK+l 


Kw^ 


(3) 


3.  Smiddy,  M.,  and  Stuart,  R.D.  (1969)  An  Analysis  of  the  Behavior  of  a Multi 

grid  Spherical  Sensor  in  a Drifting  Maxwellian  Plasma,  AFcRe-69-0013. 

4.  Vasyliunas,  V.M,  (1968)  A survey  of  low-energy  electrons  in  the  evening 

sector  of  the  magnetosphere  with  OGO  1 and  OGO  3,  J.  Geophys.  Res. 
73:2839-2884. 

5.  Formlsano,  V.,  Moreno,  G. , Palmiotto,  F. , and  Hedgecock,  P.  C.  (1973) 

Solar  wind  interaction  witti  the  earth's  magnetic  field,  1,  magnetosheath, 
J.  Geophys.  Res.  78:37 14-3730. 
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where  w is  the  most  probable  thermal  speed  and  K is  the  slope  of  the  directional 

4 

differential  flux  at  high  energies.  In  the  limit  K -►  « , this  distribution  function 
reduces  to  the  Maxwellian  case.  We  now  calculate  the  current  measured  by  a 
spherical  probe  given  by  electrons  with  a kappa  distribution. 

The  current  to  a sphere  of  radius  r is 


If 


f(v)  V 


n dSd^v 


(4) 


where  n is  an  inward  directed  unit  surface -nonnal  vector.  The  surface  element 

2 3 2 

dS  = r Sin  0 d0  d0  and  the  velocity  element  d v = v dv  SIN  e de  dX.  The  angles 

(0,^)  and  (e.X)  are  the  colatitude  and  azimuth  in  configuration  velocity  spaces 

(Figure  1)  respectively.  In  calculating  I,  the  velocity  space  integration  is  done 

1 3 

only  over  those  orbits  which  penetrate  the  grid.  * From  the  conservation  of 
energy  and  angular  momentum,  the  integration  limits  for  v and  e are 
0 < Vp  < V < 00  and  0 S e ^ e p where 

Vp  = (2  qV/m)^/^ 

and 


where  V is  the  retarding  potential  of  the  grid.  Assuming  that  the  electrons  are 
isotropically  distributed,  integration  over  configuration  and  velocity  space  angles 
gives 


I = 4^2  p2  q J 

''r 

Using  the  distribution  function  given  in  Eq.  (4)  and  allowing  for  the  grid  trans- 
mission gives 


8 


where  = l/2  m w”^  is  the  most  probable  thermal  energy.  We  note  that  in  the 
limit  K -♦  00 , Eq.  (5)  reduces  to  Eq,  (l).  An  equation  similar  to  (5)  is  derived  in 
Appendix  A for  a planar  sensor. 

As  a concrete  example  we  have  calculated  the  electron  current  to  a gridded 

_3 

spherical  electrostatic  aneilyzer  under  the  conditions  n = 100  cm  , E = 1 eV, 
and  K = 3,  5,  10.  These  conditions  could  be  realized  in  the  outer  reaches  of  the 
plasmasphere  in  the  equatorial  plane.  A short  FORTRAN  program  for  calculating 
IrCV)  with  Rq  = 3.  18  cm  and  (3  = 0.  546  is  given  in  Appendix  B.  Calculated  values 
of  log  I as  functions  of  applied  voltage  for  the  four  values  of  K are  given  in  Fig- 
ure 2.  For  simplicity,  we  have  assumed  that  the  satellite  ground  is  at  plasma 
potential.  Because  the  vehicle  potential  is  constant  during  the  retarding  portion 

g 

of  an  electron  sensor  sweep,  this  assumption  is  not  critical.  In  the  Maxwellian 
case  (K  = 00 ),  the  log  I - V curve  is  linear.  The  current  characteristic,  however, 
becomes  more  and  more  nonlinear  as  the  value  of  K decreases.  Slopes  of  the 
K < 00  curves  are  less  steep  than  in  the  Maxwellian  case.  Thus  any  attempt  to 
analyze  a K < oo  characteristic  as  though  the  plasma  were  Maxwellian,  leads  to 
an  overestimate  of  the  mean  thermal  energy  (c,  Eq.  (2)).  The  degree  to  which 
the  density  is  miscalculated  depends,  in  turn,  on  the  degree  to  which  E^  has  been 
overestimated.  It  is  thus  useful  to  develop  practical  methods  by  which  the  param- 
eters K,  E , and  n may  be  obtained  from  the  current -voltage  relation  given  in 
Eq.  (5).  We  assume  that  a In  (-I)-V  plot  has  a nonlinear  shape  indicating  that 
K < 00. 

Equation  (5)  is  of  the  form 


I„(V)  = n E*/ 

K O 


Cj  F (K) 


where  Cj  = ttR^  3 q (8a/7rm)^/^  and  F (K)  = r(K  - l)/r(K  - 1/2).  F (K) 

decreases  monotonically  from  1.  3 to  1.  0 as  K increases  from  3 to  00.  In  the  limit 

1 “K 

qV/KE^  » 1,  the  current  is  proportional  to  (qV/KE^)  . The  term  a in  the 
definition  of  Cj^  is  1.  6 X 10  erg/eV.  Thus, 


K = 1 - 


3 In  (-1) 


6,  Smiddy,  M.  , and  Stuart,  R.D.  (1969)  The  Characteristics  of  a Space -Vehicle - 
borne  Charged  Particle  Sensor,  AF(JRL-69-0519, 


y 
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Figure  2.  Semilog  Plot  of  the  Current-Voltage  Curves  to  a Spherical  Electro- 
static Analyzer  With  Rci,"  0 = 0.  546  Due  to  a Plasma  With  an  Elec- 

tron Density  of  100  cm^,  a Most  Probable  Thermal  Energy  of  1 eV  and  Spectral 
Indices  of  3,  5,  10  and  «> 


Log  (-1)  versus  log  V plots  for  the  cases  studied  above  are  given  in  Figure  3.  We 
note  that  in  the  K = 3,  5 and  10  cases  the  log  (-1)  - log  V curves  decrease  linearly 
at  large  retarding  voltages,  characteristic  of  a power  law  relationship. 

Having  determined  the  value  of  K from  the  slope  of  the  log  (-1)  - log  V curve 
at  large  values  of  qV/KE^,  we  may  now  calculate  E^.  In  the  region  V -»  0 of 
Figure  2,  the  log  (-1)  - V curve  approaches  a linear  relationship 

Lim  d In  (-1)  _ / 1 - K ^ q 

V-0  dV  ' \ K 

or  (7) 


-I  -1 

/i  - k\ 

d In  (-1) 

\ K ) ^ 

dV  J 

10 


I 


I 

I 
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Figure  3.  Log-log  Plot  of  Current -Voltage  Curves 
Under  the  Same  Conditions  Specified  for  Figure  2 


In  the  Maxwellian  approximation  (Eq.  (2))  -q  [d  In  (-I)/dV)  = kT.  Thus  if  K = 3 

the  mean  thermal  energy  is  overestimated  by  50  percent,  and  by  10  percent  if 
K = 10. 

The  electron  density  is  calculated  from  the  current  measured  when  the  sen- 
sor's grids  are  at  plasma  potential  (V  = 0) 


n = 


Cj  F (K) 


(8) 


The  error  in  the  density  calculation  is  not  so  severe  as  in  the  case  of  temperature. 
Consider  the  case  of  K = 3,  E = 1 eV  and  n = 100  given  in  Figure  2.  The  current 

-Q  ® 

at  V * 0 is  2. 42  X 10  A.  In  the  Maxwellian  approximation  F (K)  = 1 and 
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= 1.  5 eV.  is  equal  to  -1.  86  X 10  for  Rq  = 3.  18  cm  and  3 = 0.  546. 
Substitution  in  Eq.  (8)  gives  n = 106.2.  That  is,  the  density  is  overestimated  by 
6.  2 percent. 

In  Appendix  A we  show  that  for  a circular -aperture,  gridded  planar  probe 
the  current -voltage  relationship  is  the  same  as  that  given  in  Eq.  (5).  However, 
a spherical  probe  of  radius  R collects  four  times  the  current  of  a circular  planar 
sensor  of  the  same  aperture  radius,  given  the  came  plasma  parameters.  The 
above  methods  for  determining  K,  E^  and  n are,  mutatis  mutandis,  directly 
applicable. 
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Appendix  A 

Current  to  a Planar  Sensor 


The  current  across  the  gridded  aperture  of  a planar  sensor  is 
I = q //  • rt  dSd^  (Al) 

where  n = -z  is  an  inward  directed  surface  normal  unit  vector  (Figure  lb), 

f{^,  r„)  is  the  electrons'  distribution  function  at  the  aperture.  Thus  v • n = -v_„. 
o o o zo 

To  solve  equation  (Al)  we  assume  that: 

(1)  at  great  distances  from  the  aperture,  the  electrons  have  an  isotropic, 
kappa  distribution,  and 

(2 ) the  retarding  potential  distribution  is  a function  of  z,  the  distance  from 
the  z ° 0 plane,  only. 

These  assumptions  allow  us  to  approximate  the  current  to  the  sensor  as  the  pro- 
duct of  the  aperture  area  (A)  times  the  electron  flux  across  the  z « 0 plane. 

For  an  electron  to  cross  the  z > 0 plane,  at  great  distances  from  the  plane  it 
must  have  a sufficiently  high  velocity  component  in  the  z direction  to  overcome 
the  potential  barrier.  The  conditions  are 

V > Vp  ■ (2qV/m)^^^ 


15 


and 


0 S Sin  e S 


®R 


where  v and  e are  the  electron's  speed  and  velocity  space  colatitude  angle 

(Figure  lb)  at  a great  distance  from  the  detector.  From  the  conservation  of 

3 3 

energy  and  momentum,  it  is  easily  shown  that  v^^d  v^  = v^d  v.  Using  the  fact 
that  the  distribution  function  is  a constant  along  particle  trajectories  (Liouville's 
Theorem),  we  can  write 


I 


Cos  e Sin  e de 


f(v)  v^  dv 


Integration  over  velocity  and  configuration  spaces  give; 


(A2) 


2 

If  we  let  4»Rq  be  equal  to  the  area  of  the  "aperture"  of  a spherical  sensor,  then 
Eq.  (A2)  is  identical  to  Eq,  (4*)  in  the  text.  Using  the  distribution  function 


fK(v) 


_n  r(K+  1) 

w^  (xK)^'^  r{K  - 1/2) 


K+l 


gives 


Ijj(V)  ■ l/4  0 n q A c 


r(K  - 1) 
r(K  - 1/2) 


where  C = (8  E^/m)^^^,  and  0 is  the  transmission  coefficient  of  the  grid.  Note 
that  if  we  assume  that  the  aperture  is  circular  with  a radius  R^,  then  A « R^. 
Thus  for  grids  of  the  same  radius  a spherical  detector  will  collect  four  times 
more  current  than  a circular  planar  sensor. 
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Appendix  B 


C THIS  PROGRAM  13  OtSIGNEO  TO  CALCULATE  THE  CURRENT  TT  a SPHERICAL 
C ELECTROSTATI"  ANALYZER  IF  THE  AH3IENT  ELECTRONS  HAVE  A KAPPA 
C DISTRIBUTION.  ' HE  BIASES  OF  THE  GRIDS  APE  ASSUHEO  TO'flF  PETAROING 
C FOR  ELECTRON*.  I.E.  V<0. 

C THF  CURPCNT'roUATl  ON  IS  ' CURW  = *2, •OENS • Q* TRANS* ♦ SQftY 
C GAHHA(K-HE*rQRT  (K>/  (GAMMA  ( K-.  51  * ( 19>V3LT/RK*E0>  *»<M  PI 
C where"  Q=i76E-19:0UL0MBS,  TRANS  = .5<46  iTRANSHISSION  FACTOR  OF  THE 
C INJTRUHENT)  ANO  VOLT  IS  THE  APPLIED  VOLTAGE  OF  THE  GOIP.  PG  IS 
C THF  RADIUS  OR  THE  OUTER  GRID  (=.T.18CNI.  HO  IS  THF  hFAN  THERMAL  VELOCITY, 
c Eo  IS  the  mean  thermal  energy 

DIHENSTOH  VOLT'Ci  3»,A1)D(  13)  ,0  NMCiX),  CURVE  13 ) 

_1  REAOtS.in)  OENS.EO.K 
10  FORMATEZriJ. 1,151 
PI  = 3.  ILI'593 
"pi  12=SORT< 31) 

R<=1C 

"KHOik-1 

EO12=S0R’<E0) 

C CALCULATE  GA Hl= G A MH A <K- 1 ) = (K-2)) 

P.N=  <-?. 

TF(RN)  Z,7.A 

7 GA-lrl. 

GO  TO  <j 

8 GAH1=RN 

9 PN=RN-1. 

IF  (RN)  7,8.5 
5 GAN1=GAH1»RN 
GO  TO  9 
6''C0NtlNUE 
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IBIS  PACE  IS  BEST  QUALITY  PBACTICABL^ 

JBQM  OQPY  piaaiisaaD  TO  BQQ  ^ 


C OLCULATE  GAN?  = GftHHA<K-l/2)  * SOR^LPIl.*  < 2<"3U  • 1 • 

L = 3 

NUM=2»<-1 

18  L»L»2 

NUM1=2»<-L 

IF<NUMl-n  16,16,1/ 

17  NUHaMUH^NU*^! 

GO  TO  18 

Ife  GAM2sNUi»Pri2/ (2 .** ( K-l»  » 

FACT=GAH1*SJPT (RK)/GAM2 
WorTE(6.?9»0ENS,C0tRK 
99  FORMAHSFll.l) 

^0.50J<  = 1,5 

no  A9  J=l,13 

GO  TO  (21,22,23,24,291  K 

21  VOLTC J»=(J-1»*0. 001 

GO^  TO  30 

23  tfOLTC J»  =■( J-1 
GOjro  30 

22  VOLT(j»=(J-ll*0.01 
___6_0.T0  30 

24  VnLT(J)=(J-l»*l. 

GO  TO  3 0 

29  VOLT(J»=(J-l J*10.0 

30  AOOCJIsVrLTC Jl /(PK*EO» 

P^4MtJ^  = (l.♦AOP(J^»**K^10 

CURV( Jl=l. 9aE-10*DENS*SQRT (E 0) ♦ FACT/ONK< J» 

WRITE(6j48) #OLT( J»,CURVIJ» 

48  FORMAT«2F20.9» 

_4a.  CONTINUE 
90  CONTINUE 
_ GO  TO  1 
END 
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